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Using purely phenomenological approach we show that the peaking structures observed in the 
ip(2S)n + and Xci 7r+ mass spectra in B — > %p(2S) (xci)tt + K decays can be result of (DD ( -*^) + — > 
(cc) Ics n + rescattering in the decays B D~j(— > D^K)D^*\ In particular, the position of the 
peak in the chain B — > D S (2S)~ D + — > K~ D D + — > K ~ip( 2 S)tv + coincides well with the measured 
Z(4430) + mass, assuming the mass of D S (2S)~ (the first radial excitation of DJ) to be 2610 MeV/c 2 . 
The widths of the Z(4430) + peak is also well reproduced in this approach independent on the 
width of D S (2S)~ . Although the decay B — > D S (2S)~ D + has not been observed so far and even 
^-j, D S (2S)~ -meson is not discovered yet, this decay is expected to be large, and the mass of D S (2S)~ is 

predicted in the range (2600-2650) MeV/c 2 . The broad bump in Xci 7I " + spectrum can be attributed 
to the B — > D S (2S)*~ D + — ► K~ D°D + decay observed with a large branching fraction followed by 
rescattering D°D + —> Xci 7r+ - 
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. The charmonium-like charged Z + states, seen by Belle in B — > ip(2S)n + K and B — > Xci^ + K decays remained 
puzzles for last few years. Belle [l[ observed the first Z(4430) + state as a sharp peak in ip{2S)Ti + mass spectrum near 
O !■ M(ip(2S)ir + ) — 4430MeV/c 2 with statistical significance of more than 6 a. The main background around the peak 
region is due to B -> tp(2S)K*(892) and B -> ip(2S)K% (1430) decays. In the Belle analysis the two if*W states were 
vetoed to suppress this background. It was noted that interference between different partial waves in the 7rif-system 
can produce peaks, that are reflections of the K*(*' polarization. However, these effects should also produce additional 
y—i | sharp structures nearby in M(if)(2S)ir + ), which are not observed by Belle. More detailed analysis was performed by 
Belle a year latter to prove quantitatively the absence of fake peaks by the Dalitz fit over all signal events including 
regions. In this study Belle confirmed the Z(4430) + observation and found its parameters consistent with those 
measured in the first paper. Although BaBar Q has reported no evidence of Z(4430) + in their B — > i/;(2S)ir + K 
analysis, the non-uniform structures are presented in their spectrum as well, and both Belle and BaBar spectrum are 
consistent with each other. 



Two broader peaks (Zi(4050)+ and Z 2 (4250)+) were found by Belle 0] in the Xci 7r+ mass spectrum in the analysis 
of B — > Xci^ + K° decay. A Dalitz fit with a single resonance in the Z + channel is favored over a fit with only 
I ■ K *(*> -resonances and no Z + -fit by more than 10 a. Moreover, a fit with two Z + resonances is favored over the fit 
• • ' with only one resonance by 5.7 a. 

If the observed Z + peaks are real states, they would necessarily be exotic (non-conventional qq) mesons, as their 
minimal substructure consists of four quarks. Many attempts of explanation of Z + 's follow these observations including 
molecular 0, Q , tetraquark 0| , hadrocharmonium states [1] or cusp effects Q . In this Letter we demonstrate that 
■ the observed peaks can be explained by the effect of rescattering in the decay chain 

B D'jDW followed by ->• D&K (1) 

of D(*)£>W -pair into charmonium+7r + , where one is directly produced in B-decay, while is from the 

intermediate D~j resonance (Fig. [T|) . 

We assume that the B decay dynamics can be factorized from the rescattering process. Under this assumption 
the mass of the charmonium+7r + combination is equal to those of D^*'D^*' . The mass spectrum of (cc) rcs tt + system 
produced via rescattering consists of peaking structure(s) reflecting the D~ } polarization in B-decay or D^D^ 
polarization in the formation of (cc) rcs 7r + . Here we denote the D^D^ = (cc) res 7r + as "Z" and calculate "Z" mass 
ignoring the subsequent decay 11 Z" — > (cc) Tes ir + . However, angular momentum/parity conservation impose restrictions 
on "Z" production depending on the final (cc) ros 7r + state that should be taken into account. In particular, Tp(2S)ir + 
and Xci 7r+ systems have different spin-parity, thus different decay chains of the type ((TJ should be proposed for the 
explanation of the peaks in their spectra. We also assume that rescattering of D^D^ in S*-wave dominates. 

There are many decay chains ([1]) that can provide the required conditions; all of them should be considered taking 
into account interference between them. We are now limiting ourselves by searching for the dominant contribution 
that can roughly reproduce the features of the observed (cc) ros ir + spectra. We note that orbital D~ excitations are 
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FIG. 1: Rescattering diagrams resulting in B — ¥ (cc) rcs k + K final state. 



TABLE I: The angular parts of the B -> D~jD w ->■ {D M D (t) )K ->■ "Z"if chain matrix element. 
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hardly suitable for our explanation: j — 1/2 states are below D^K threshold, while B-decays into j = 3/2 states are 
strongly suppressed (£> < 10~ 4 [13]). Radial Dj" excitations could be better candidates: the decay B — > D s (2700)~ D 
was observed by Belle with relatively large branching fraction B(B -)• D s (2700)-D) x B(D S (2700) _ ->■ DiC) ~ 
10~ 3 [ll|. The quantum numbers of the Z? s (2700)~ have been measured to be J p = 1~ and this state is likely to 
be a radial excitation of D* s " (D S {2S)*~). If this is really true, the pseudo scalar state, D S (2S)~, have a mass of 
(2600 - 2650) GeV/c 2 (expected 2S 1 - 2S 3 splitting is (60 - 100) MeV/c 2 [H), and should decay predominantly into 
D*K. It is also expected that two body -B-decays into D S (2S)~ are not suppressed. 

We are first looking for an explanation of the Z(4430) + peak. The rescattering of {DD) + — > i[i(2S)ir + is forbidden 
by parity/angular momentum conservation, while the diagram (DD*) + — > ^{2S)tt + is allowed for both initial and 
final systems in the S'-wave. The two B — > D + D*°K and — > D°D* + K decay modes that include all intermediate 
states (as expected mostly D~j) has a relatively large branching fraction of ~ 10 -2 [l(|. (DD*) + system should form 
a pseudo-state with the spin equal to one and positive parity ( J p = 1 + ), thus the parity of the intermediate D~jD^ 
combination should be also positive. The decay chain B — > D S (2S)~~D + followed by D S (2S)~ — > D*°K~ matches 
this parity constraint. Another allowed decay chain that can result in i(j(2S)tt + is B — > D S (2S)*~ (— > D°K)D* + in 
case of S (-D)-wave between D S (2S)*~ and D* + . We calculate the matrix elements of these decays in the helicity 
formalism 
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where Abw is a Breit-Wigner function for a corresponding D s j-resonance; D{6, \d*) is it's angular part depending 
on the D~j decay angle, 9, and helicity of D*, (In case of the second decay chain the angular term depends on 

the D S (2S)*~ helicity; however here we use equality \d*(2S) = Ajj*). The next term, -Dd* (0', \d* ), is responsible for 
the rotation of the D* spin from the D~j (or B in case of the second decay chain) to "Z" rest frame. We note that in 
the 11 Z" rest frame the D* helicity is fixed to zero. Finally, Dz{9") provides formation of "Z" from helicity-0 D* and 
D. We list explicitly angular contributions to the matrix elements in the upper lines of Table U for two considered 
decays: B — > D S (2S)~D and B — > D S (2S)*~ D* in terms of (small) Wigner functions. In the later decay we assume 
that S (over D) wave dominates. 

We calculate the M("Z") spectra using Monte Carlo simulation. In this calculations we assume the mass of 
A i (2S')- to be 2610MeV/c 2 and its widths to be 60MeV; the D S (2S)*- parameters are fixed to PDG values 
The obtained spectra are presented for B —> D S (2S)~ D and B — > D S (2S)*~ D* decays in Fig.[2]a) and b), respectively. 
The solid lines show the Z(4430) + peak position; the dashed lines show ±T Z + window (as determined in Ref. [2||). 
The main peak position and width in Fig. [2] a) are close to the experimentally measured values for Z(4430) + . The 
contribution from B — > D S (2S)*~ D* (Fig. [5] b) may be responsible for another broad excess of signal events over 
combinatorial background in the region above 4.1GeV/c 2 , that can be seen in both Belle and BaBar data. For 
comparison the experimental spectra seen at Belle and BaBar are shown in Fig. [5]c) and d), respectively. We note 
that the instrumental reconstruction efficiency in the region of (4.6 — 4.75) GeV/c 2 drops very sharply, as this mass 
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FIG. 2: The M("Z") spectra in the decays a) B -> D S (2S)~D; b) S -> D S (2S')*"D*. The ip(2S)-K + spectra in the selected 
B V>(25)7r+A' events in the c) Belle and d) BaBar data. 



region corresponds to a low center-of-mass momentum of kaon. Thus experimental bias can hide the accompanying 
high mass structure which can be seen in our Monte Carlo spectra. 

For the explanation of the broad bump in Xci^ + spectra around 4 — 4.4 GeV/c 2 region another chain of the type (p} 
with negative parity of the final state has to be proposed. At least one P-wave is required in the rescattering 
(_D(*)£)(*))+ — s> Xci 7r+ to provide parity conservation. The simplest one (DD)t, — > (xd^ + )p can be provided by 
the known decay B — > D S (2S)*~ D. The Monte Carlo (DD) + mass spectrum from this decay is shown in Fig. [3] 
a). The angular part was generated according to the elements listed in Table fl] For comparisons the Belle Xci^ + 
spectrum from B — > Xci^ + K decays (from Ref. [1]) is shown in Fig. [3]b). We remind that the experimental efficiency 
around high mass bump is significantly lower, which can explain much smaller peak in the experimental spectrum at 
(4.6 - 4.8) GeV/c 2 . 

In summary, we show that the {DD*) + — > ip(2S)ir + rescattering in the decay chain B — > D S (2S)~ D, D S (2S)~ —> 
D*K can explain appearance of a peaking structure in ij)(2S)ir + mass spectrum in B — > tp(2S)7r + K decays around 
4.43GeV/c 2 . The broad structure in Xci 7r+ mass spectrum around 4 — 4.4GeV/c 2 can be also explained within 
the similar approach using another decay chain B — > D S (2S)*~ D, D S (2S)*~ — > DK with subsequent rescattering 
(DD) + — > Xci7r + - We note, that if our explanation is valid for charged charmonium-like peaks observed in £?-decays, 
the similar mechanism should also reveal itself in other processes. For example, large e + e~ — > ip(2S)ir + Tr~ cross 
section can be also explained by rescatteringfrom the e + e~ — > DD*tt process followed by DD* — > i[i(2S)ir. The later 
cross section have been measured by Belle [13] , though with large uncertainty. The peaks in the former cross section 
(considered now as the F(4350) and F(4660) states) may be related to the interference of ^-resonances and many 
intermediate states (e.g. DiD, D2D, D2D*, etc.) in their decays with different D* helicities resulting in the DD*ir 
final state. In particular, large e + e~ — > -0(4415) — > DD2 cross section [lj] can serve as a hint of the vicinity of the 
y(4350) and F(4660) states to the '0(4415) resonance. However, the quantitative calculations are still beyond our 
capabilities mainly because of lack of experimental data. 
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